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Abstract 

For static reductions of isotropic and anisotropic Magnetohydrodynamics plasma equilibrium 
models, a complete classification of admitted point symmetries and conservation laws up to 
first order is presented, ft is shown that the symmetry algebra for the isotropic equations is 
finite-dimensional, whereas anisotropic equations admit infinite symmetries depending on a free 
i-C . function defined on the set of magnetic surfaces. A direct transformation is established between 

isotropic and anisotropic equations, which provides an efficient way of constructing new exact 
, anisotropic solutions. In particular, axially and helically symmetric anisotropic plasma equilibria 

' arise from classical Grad-Shafranov and JFKO equations. 
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1 Introduction. 



Systems of isotropic Magnetohydrodynamics (MHD) and anisotropic Chew-Goldberger-Low (CGL) 
plasma equations, in particular, their equilibrium reductions, are used for description of plasmas in 
controlled thermonuclear fusion research, geophysics and astrophysics (Earth magnetosphere, star 
^ ■ formation, solar activity), and laboratory and industrial applications [1—3]. 

MHD and CGL systems, as well as their equilibrium versions, are essentially nonlinear systems 
of partial differential equations in 3D space. Knowledge of physically meaningful exact solutions 
and analytical properties of these systems (such as symmetries, conservation laws, stability criteria, 
etc.) is important for understanding the core properties of the underlying physical phenomena, for 
modelling, and for the development of appropriate numerical methods. 

Common ways of finding exact analytical solutions to such systems include reduction by a symme- 
try group (similarity solutions), the use of symmetry transformations to generate new solutions from 
known ones, and the use of mappings from solutions of other equations. The first approach applied 
to axially and helically symmetric static MHD configurations has yielded the well-known Grad- 
Shafranov [4-6] and JFKO [7] equations, and led to several classes of exact solutions (e.g. [8-11]). 
A different approach that makes use of equilibrium solution topology (general existence of 2D 
magnetic surfaces) was used in [12]. 
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A symmetry of a system of PDEs is any transformation of its solution manifold into itself. A sym- 
metry thus maps any solution to another solution of the same system. Several types of symmetries, 
such as continuous (point, contact, higher-order, nonlocal) Lie groups of symmetries and discrete 
symmetries, can be obtained algorithmically (e.g. [13-17]). In particular, using Lie's algorithm for 
solving symmetry determining equations, one can discover one-parameter, multi-parameter, and 
infinite-dimensional symmetry groups. Symmetries are used as transformations that yield new so- 
lutions and new conservation laws of differential equations from known ones, and also for finding 
particular symmetry-invariant solutions. Knowledge of symmetries is essential to answer the ques- 
tion about the possibility of mapping a given PDE system into a target PDE system or a class of 
PDE systems. 

An important complement to the full symmetry structure of a PDE system is knowledge of its 
conservation law structure. Local conservation laws contain important information about physical 
properties of a model under consideration, and provide conserved norms used in analysis of solutions 
and also in development of numerical methods. Conservation laws can be found algorithmically by 
a direct construction method in terms of multipliers that satisfy determining equations related to 
the adjoint of the symmetry ones [18,19], without the need for any Lagrangian. In particular, this 
method allows one to by-pass all limitations of Noether's theorem. 

The paper is organized as follows. In Section [2j we describe the static plasma models under 
consideration, as well as their basic properties, and state the transformation that relates the two 
models. In Section [3l we classify and compare point symmetries of MHD and CGL static plasma 
equilibrium systems. In particular, we demonstrate that the static CGL system admits an infinite- 
dimensional symmetry group (which appears to be related to the infinite-dimensional symmetry 
group of dynamic CGL equilibrium system [20]). We derive infinitesimal and global representations 
of the admitted symmetry groups, and discuss physical properties and group structure of the 
infinite symmetries that arise for the static CGL system. In Section we discuss applications 
of the infinite symmetry group to construction of exact static anisotropic plasma equilibria. In 
particular, we show that axial and helical static anisotropic (CGL) equilibria arise from solutions 
to conventional Grad-Shafranov and JFKO equations. An explicit example of an exact solution 
describing an anisotropic plasma vortex is presented. Finally, in Section [5l we complete the analysis 
by classifying all conservation laws admitted by static MHD and CGL systems with multipliers 
linear in first-order partial derivatives, and discuss their physical meaning. From the comparison 
of symmetry and conservation law classifications of these two static plasma equilibrium systems, 
we establish a direct transformation (see Theorem [TJ) between the plasma models (including the 
equivalence of solution sets.) In Section [6l we summarize the results presented in this paper in a 
larger context of dynamical MHD models. 

Symbolic software packages GeM for Maple [21] and Crack/LiePDE/ConLaw for REDUCE [22] were 
used for all symmetry and conservation law computations. 

2 Static plasma equilibrium models 

Equilibrium (time-independent) plasma models with and without flow are used in many physical 
applications. In particular, for the purpose of analysis, static plasma equilibrium systems are often 
considered. On one hand, static equilibrium equations are much simpler than dynamic ones, and 
yield to analytical techniques more easily; on the other hand, they are still nonlinear 3D models, 
which inherit many properties from full plasma equilibrium models. Static plasma equilibria, even 
in a simplified force-free (constant-pressure) setting, are used in astrophysical modelling. 

The static MHD equilibrium equations are obtained directly as a static time-independent reduc- 
tion of the full system of MHD equations (e.g. [1]). They have the form 



curl B x B = grad P, div B = 0. 
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Here B is the vector of the magnetic field induction, and P is plasma pressure. The electric current 
density is given by J = curlB. 

In all static isotropic plasma equilibria (|2.ip . the magnetic field B is tangent to 2-dimensional 
magnetic surfaces *f?(x,y,z) = const that span the plasma domain: B • grad^ = 0. The plasma 
pressure P = P(&) is constant on magnetic surfaces. In a compact domain, magnetic surfaces are 
generally tori [23]. When magnetic field lines are closed loops or go to infinity, magnetic surfaces 
may not be uniquely defined, and one may specify ^(x, y, z) = const on each magnetic field line. The 
only possible case when magnetic surfaces do not exist is Beltrami- type configurations curlB = aB, 
a = const. 

We note the well-known equivalence between the static plasma equilibrium system (|2.ip and the 
time-independent Euler's equations of inviscid fluid motion 

(v • grad) v = — grad tt, div v = 0, (2-2) 

IBI 2 

with velocity v = B, pressure tt = Pq — P — and constant density p = 1. 

The static time-independent equilibrium of anisotropic plasmas is a similar reduction of the CGL 
equations [2,20,24] 

|B| 2 

(1 — r) curl B x B = grad p± + r grad— h B(B • grad r), div B = 0. (2.3) 

Here r is the anisotropy factor, and pressure P is a symmetric tensor with two independent param- 
eters p\\,p±: 

r-a^ft. P-Bj + aZftBaB, (2 .4) 

where I = diag(l, 1,1) is the identity tensor in M 3 . 

Equations (|2.3p describe equilibria of strongly magnetized or rarified plasmas; py denotes pressure 
along the strong magnetic field B, and p± is pressure in the transverse direction. Unlike static 
isotropic plasma equilibria, anisotropic plasmas described by (|2.3p in general do not possess 2D 
magnetic surfaces. 

Remark 1. Static plasma equilibrium systems (|2.ip and (12. 3D arise from Boltzmann and Maxwell 
equations under essentially different isotropy assumptions [2,25,26]. In particular, the isotropic 
model ()2.ip was derived from Boltzmann equation using expansion in powers of mean free path, 
while for the anisotropic model the expansion in powers of ion Larmor radius was used, which 
implies py / pj_. However it is easy to see that for r = equations (12. 1|) and (12. 3|) coincide. 

Closure of the anisotropic equilibrium system (|2.3p requires an equation of state. Note that 
dynamic equations of state, such as "double-adiabatic" equations 



d ( P± \ d ( P|||B| 



dt Vp|B|y dt \ p 
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suggested in the original CGL paper [2], cannot be used since they vanish identically for all static 
equilibria. In this paper we use the equation of state [20] 

B-gradr = 0, (2.5) 

which implies that the anisotropy factor r is constant on magnetic surfaces (more generally, on 
magnetic field lines). Below we show that the relation (|2.5h can be rewritten in a form (12. 8p that is 
similar to the equation of state used in numerical modeling of anisotropic magnetosheath plasma 
in the CGL approximation [27] (see also [12]). 
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In Sections [3] and [5j we classify and compare point symmetries and conservation laws admitted by 
the isotropic (MHD) static plasma equilibrium model (|2.ip and the anisotropic (CGL) model (|2.3p 
with the equation of state (|2.5p • The point symmetry classification of the anisotropic model is found 
to differ from that of the isotropic model by only one symmetry generator (|3.1ip depending on 
an arbitrary function defined on magnetic surfaces, which describes infinite symmetries admitted 
by the anisotropic system (|2.3p . (|2.5p . 

Since any solution (B,P) of static isotropic (MHD) equilibrium is automatically a solution 
(B,pj_,r) = (B,P, 0) of the static anisotropic (CGL) system (|2.3p . (|2.5p . the infinite symme- 
tries Yqq map each MHD equilibrium solution into a continuum of CGL equilibria. This suggests 
that a direct relation may exist between the two static plasma equilibrium systems. Comparison 
of multipliers and fluxes of conservation laws admitted by the two plasma equilibrium systems (see 
Section [5]) indicates the specific form of such a transformation: the isotropic equilibrium magnetic 
field B should be proportional to y/1 — rB of the anisotropic model. This leads to the following 
important theorem. 

Theorem 1. The static anisotropic (CGL) equilibrium system 112. 3]) with equation of state \2.5\) 
can be written in the form 

curl (VI - rB) x (Vl - rB) = gradp, (2.6) 



div(Vl - rB) = 0, (2.7) 
B • gradr = 0, 

where p = p± + ^rB 2 = \{p\\ + p±) is the mean pressure of anisotropic plasma configuration. 
This form of the static anisotropic (CGL) equilibrium system A2.3]) . 112.5]) is equivalent to the static 
isotropic (MHD) equilibrium system h2.1]) with pressure p, magnetic field y/l — rB, for any smooth 
function r constant on magnetic surfaces / magnetic field lines. 

The proof proceeds by expanding (|2.6p . (|2.7p using vector calculus identities. 

Corollary 1. The 3D domain of every static anisotropic (CGL) equilibrium plasma configuration 
satisfying (12. 3|) . (|2.5p is spanned by 2D magnetic surfaces 

y, z) = const 

(with the only possible exception being Beltrami-type configurations curl y/1 — rB = ct\/l — rB, 
a = const. ) In particular, the mean pressure p = p(^) and the anisotropy factor r = t(*$>) are 
constant on magnetic surfaces. 

If magnetic surfaces are not uniquely defined (magnetic field lines are closed loops or go to infin- 
ity), *If(x,y,z) can have a different constant value on every magnetic field line. 

This corollary follows from the equivalence of the equations (12.6P to the static MHD system (12. ip . 
for which this result is well-known. 

We also note that since r = r(^), the equation of state (|2.5p can be rewritten as 

P±/p\\ = 1" (rW|B| 2 /P||) • (2.8) 
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3 Point symmetries of plasma equilibrium models 

We will classify all point symmetries admitted by static isotropic (MHD) and anisotropic (CGL) 
systems. In particular, we seek Lie groups of point symmetries admitted by the system (|2.ip 



x' = x + e£(x,B,P) + 0(e 2 ), 

B' = B + eq(x, B, P) + 0(e 2 ), (3.1) 
P> = P + ea(x,B,P) + 0(e 2 ) 

with corresponding infinitesimal generator 

where the hook denotes summation over vector components. Symmetry components (£, rj, a) satisfy 
determining equations given by the invariance of the solution set (B(x),P(x)) under X: 

curl B x (t] — • grad)B) + curl (rj — (£ ■ grad)B) x B = grad (a — (£ ■ grad)P), , , 
div (ry - (£ • grad)B) = [ } 

holding for all static equilibria satisfying (12.11) . 
Similarly, a symmetry generator of static anisotropic equilibrium system (|2.3p . (|2.5p has the form 



+ ^ 9 d d 

whose components (^, 77, cjj_, v) are found from corresponding determining equations that state the 
invariance of the solution set of the equations (|2.3p , (|2.5p under the action of Y (j3.4j) . 



3.1 Point symmetries of static MHD and CGL plasma equilibrium systems 

The following theorems describe all Lie point symmetries of static isotropic and anisotropic plasma 
equilibrium PDE systems. The proofs are computational and consist of applying Lie's algorithm 
[13, 14] for solving the determining equations (|3.3|) and their counterpart for the anisotropic case. 

Theorem 2 (Point symmetries of static MHD equations). The point symmetries admitted by 
the static isotropic plasma equilibrium system (|2.ip form a nine- dimensional Lie algebra with the 
following generators: 

• Killing symmetries (translations and rotations) 

Xk = C-'tt" + (B • grad)£j— — ■; (3.5) 
ox oa 

• Scalings and dilations 

• Pressure shifts 

* = I- (3-7) 
Here ( = a + bxxisa Euclidean Killing vector; a,bgK 3 are constant vectors. 
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Theorem 3 (Point symmetries of static CGL equations). The static anisotropic plasma equilibrium 
system (|2.3p with equation of state (|2.5p admits an infinite- dimensional symmetry algebra spanned 
by the following generators: 

• Killing symmetries (translations and rotations) 

yK = Cj ^ + (B " grad)Cj M 5 (3 ' 8) 

• Scalings and dilations 

d d d 

y s = Bj— + 2 P± — , Fd = xj- ; (3.9) 
aB op± ox 

• Pressure shifts 

Y P = (3.10) 
d P ± 

• Infinite- dimensional transformations 

^ = /B J A_ / |B|^ + 2/(l- T )| : . (3.11) 

Here ( = a + bxxiso Euclidean Killing vector (a., b E M 3 are arbitrary constant vectors), and 

IBI 2 

/ = / (p, t) is an arbitrary smooth function (p = p± + T 1 —^- is the mean pressure of anisotropic 
plasma configuration). 

The point symmetry classifications of the static MHD and CGL equilibrium systems presented 
in Theorems [2] and [3] coincide except for the infinite symmetries (|3.1ip of the CGL system. 
The existence of these infinite symmetries underlies the equivalence of the static MHD and CGL 
equilibrium systems stated in Theorem [TJ Corollary [1] of Theorem [H in turn, leads to the following 
clarification of the form of arbitrary function in the infinite symmetry generator Y^. 

Corollary 2. The arbitrary function f = f (j>, t) in infinite symmetries (13. lip of CGL equilibrium 
equations can be expressed as a function f = of a single variable ^ = ^(x) which enumerates 

magnetic surfaces (or, in general, magnetic field lines) of the original plasma configuration. 

3.2 Finite form and other properties of symmetry transformations 

Using a standard reconstruction formula (Lie's First Theorem) [13,14], we find the global Lie groups 
of point transformation groups corresponding to symmetry generators (|3.5p - (|3.1ip . 

1. Translations (Xk,Yk)- 

x' = x + a (a = (oi, 02,03) G BL 3 ). 

2. Rotations (Xk,Yk), parameterized, for example, by Euler angles <j),6,tp: 

x' = A 3 A 2 A^; B'(x') = A 3 A 2 A 1 B{ J c), 

cos (j) sin <^> \ / 1 0\ / cos ^ sm ^ 

v4i = I — sin <p cos 4> ; A 2 = cos 9 sin 6 \ ; A3 = l — sin tp cos ip 
1/ \ -sinfl cos<9 / V 1 
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3. Scalings (X S ,Y S ): 



B' = a 4 B, p' = 2a 4 V (o 4 el). 
(V = P or p± for MHD and CGL respectively.) 

4. Dilations (X D ,Y D ): 

x' = C15X (05 e M). 

5. Pressure shifts (Xp, Yp): 

p' = p + a 6 (a 6 el). 
(V = P or p ± for MHD and CGL respectively.) 

6. Infinite-dimensional transformations (l^o, anisotropic equilibria only): 

B' = M(*)B, t' = 1 - (1 - r)M" 2 (^), 

A =P±+ B2 " 2 (B ° 2 , P|| = Pi + (B') 2 (l - (l - ^-^) M- 2 (*)) , (3.12) 

where ^ is a function enumerating magnetic surfaces (more generally, magnetic field lines) of the 
original plasma configuration (B,pj_, py), and the arbitrary function M(ty) is related to the function 
f(p, r)=/(#) in HSU}. 

It is easy to see why no infinite transformations similar to Y^ are admitted by the static isotropic 
plasma equilibrium system (|2.ip . in spite of the equivalence (Theorem [1]) between the isotropic and 
anisotropic systems. Indeed, from (13. 12ft it directly follows that for anisotropic equilibria, quantities 
y/1 — rB and p± + ^tB 2 are invariant with respect to YJxj. Hence, according to the relations (|2.6|) . 
(|2.7p . these infinite symmetries correspond to the identity transformation for isotropic plasma 
equilibria. 

The structure and properties of infinite transformations f|3. 12j) and resulting families of solutions 
are discussed below. 

3.3 Properties of the infinite symmetries (13.121) of the static CGL system 

The infinite symmetries (|3.1ip . (|3.12p of the static CGL system constitute a subgroup of the infinite- 
dimensional symmetry group G of dynamic CGL equilibrium system found in [20], and share many 
of their properties. 

(i) Structure of the arbitrary function. The transformations (|3.12p depend on the topology 
of the solution they are applied to, namely, on the set of magnetic field lines = const. For the 
following topologies of the initial solution, the domain of the function ^ = ^(x) is evident: 

1. The lines of magnetic field B are closed loops or go to infinity. Then the function ^(x) is 
constant on each magnetic field line. 

2. The magnetic field lines are dense on 2D magnetic surfaces spanning the plasma domain T>. 
Then the function ^(x) has a constant value in on each magnetic surface, and is generally 
a function defined on a cellular complex (a combination of 1-dimensional and 2-dimensional 
sets) determined by the topology of the initial solution {B,pj_,p||}. 

3. The magnetic field lines are dense in some 3D domain D. Then the function ^(x) is constant 
in V. 
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(ii) Topology, boundary conditions, and physical properties. From (|3,12p it is evident that 
B' || B, thus magnetic field lines of the original plasma configuration {B,p_|_, py} are retained by 
the transformed solutions {B',pj_,pj|}. Therefore usual plasma equilibrium boundary conditions of 
the type n • B|#p = (n is a normal to the boundary &D of the plasma domain) are preserved. 

If the free function M(^f) is separated from zero, the transformed solutions retain the boundedness 
of the original solution; the same is true about the magnetic energy B 2 /2. For models in infinite 
domains, the free function must be chosen so that new solutions {B,pj_, py} have proper asymptotic 
behaviour at |x| — ► oo. 

(iii) Stability of new solutions. No general stability criterion is available for MHD or CGL equi- 
libria. However, several explicit instability criteria are known. In particular, under the assumption 
of double-adiabatic behaviour of plasma [2] the condition for the fire-hose instability is [28] 

p {l - P± >B 2 , (3.13) 

or, equivalently, r > 1. According to (|3.12jl . we have 

1-ti = (l-r)/M 2 (*), 

hence the infinite transformations A3. 12H do not change fire-hose stability /instability of the original 
plasma configuration. 
The mirror instability [28] occurs when 

It can be shown that for every initial configuration {B,pj_,p||}, there exists a nonempty range of 
values which M(\P) may take so that the mirror instability does not occur. The proof is parallel to 
that in [20]. 

(iv) Symmetry group structure. We consider the set Gc of all transformations f|3. 12 j) with 
smooth M(\P). Each such transformation is uniquely defined by a pair {a, H(*S>)}: 

M (¥) = a exp (#(¥)), a = ±1. 

The composition of two transformations (a, H and (/?, K{^)) is equivalent to a commutative 
group multiplication 

and the inverse (a, H)^ 1 = (a, —H). The group unity is e = (1, 0). Thus Gc is an abelian group 

G c = A* Z 2 (3.15) 

with two connected components; ^4^, is the additive belian group of smooth functions in M 3 that 
are constant on magnetic field lines of a given static CGL configuration. 



4 Construction of anisotropic plasma equilibria 

In this section, we show how the infinite group of transformations (|3.12p is used to construct families 
of anisotropic (CGL) plasma equilibria from a single known exact static CGL or MHD solution. 
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4.1 Construction of general 3D anisotropic plasma equilibria 

Theorem 4 (Construction of anisotropic plasma equilibria). For any given solution {B,pj_,p||} 
of the static anisotropic (CGL) plasma equilibrium equations (|2,3p . (|2.5p . or any given solution 
{B,P} of the static isotropic (MHD) plasma equilibrium equations (|2.1|) . there exists an infinite 
family of anisotropic (CGL) plasma equilibrium solutions given by (|3.12p depending on an arbitrary 
function of one variable. The infinite family of solutions has the same set of magnetic field lines 
as the original solution. 

This theorem follows directly from Theorem Q] and infinite transformations (|3.12|) . 

In particular, for any given solution {B, P} of the static isotropic (MHD) plasma equilibrium 
equations (|2.1|) . from (|3.12|) we see that the corresponding infinite family of anisotropic (CGL) 
plasma equilibrium solutions {B',pY)Pj|} is given by 

B' = M(*) B, t' = 1 - M~ 2 (*), 

p' x = P 1 + P + hi 2 (1 - M 2 (M>)) , p > { =P 1+ P- I B 2 (1 - M 2 (*)) , (4.1) 

where P\ is an arbitrary constant, M(^) is an arbitrary smooth function, and ^ = ^(x) enumerates 
magnetic surfaces (or, in general, magnetic field lines) of the original isotropic plasma configuration. 



4.2 Axially and helically symmetric anisotropic plasma equilibria 



(i) Grad-Shafranov equation. Bragg and Hawthorne [4] in 1950, and Grad and Rubin [5] and 
Shafranov [6] in 1958, have shown that the static isotropic (MHD) system f)2. 1 1) with axial symmetry 
(independent of the polar angle (ft) is equivalent to one scalar equation, called Bragg-Hawthorne or 
Grad-Shafranov equation (GS): 

^rr — — + + /(*)/'(*) = -r 2 P'(*), (4.2) 
r 

Here (r,4>,z) are cylindrical coordinates, ^(r, z) the unknown flux function (the function enumer- 
ating magnetic surfaces), and pressure P(^f) and function I{^) (related to poloidal magnetic field) 
is an arbitrary function. Primes denote derivatives. 
The expression for the magnetic field B is 

^ z I(^) $ r , s 

B = -^e r + -^e^ - -^e z . 4.3 



(ii) JFKO equation. Another reduction of static isotropic plasma equilibrium equations (|2.ip 
describes helically symmetric plasma equilibrium configurations, i.e. configurations invariant with 
respect to the helical transformations 



z — ► z + 7/1, (p — ► ip + h, 



r. 



Equations (|2.ip in helical symmetry also reduce to one equation (Johnson- Frieman-Kruskal-Oberman, 
or JFKO equation [7]) for one unknown function ^(r, u), where u = z — 7^ is the helical coordinate, 
and 7 = const 7^ 0: 



j.2 f 



f2 _|_ ^y2 



_|_ /-y2 



2 7 I(*) 
(r 2 + 7 2 ) 2 



-/iP'(tt). 



(4.4) 
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At 7 = 0, (|4.4p reduces to the GS equation (|4,2p . Again magnetic surfaces are given by \&(r, z) = 
const, and P(^f) is plasma pressure. The helically symmetric magnetic field is given by 

rl(^f) 7/(^) - r^ r 

B = -Av + 5^ + 5.6,, 2?, = ^ ] \ B z = 7 ^ r . (4.5) 



(iii) Axially and helically symmetric anisotropic (CGL) plasma equilibria. Though 
an analogue of the GS equation exists for anisotropic (CGL) plasmas [29], this equation is so 
complicated compared to the original GS equation, that using it for finding exact solutions is 
impractical. 

A practical way of obtaining families of axially and helically symmetric anisotropic plasma config- 
urations (12. 3p is through the application of transformations (14. lj) to known exact solutions of con- 
ventional GS or JFKO equations (many of such solutions are known in literature, see e.g. [8,9,11]). 
For a given axially or helically symmetric isotropic (MHD) plasma equilibrium (B, P) in a domain 
P, transformations (|4,ip yield a family of anisotropic (CGL) equilibria. This family of solutions in- 
volves an arbitrary function M(\&) defined, in general, on a celluar complex, as follows (see Section 
ES part (i)): 

• In parts of the plasma domain T> where magnetic field lines are dense on closed 2D magnetic 
surfaces, M(^) is a function of one variable \& enumerating magnetic surfaces; 

• In parts of T> where magnetic field lines of the given MHD equilibrium are closed loops or 
go to infinity, values of M(^) can be chosen independently on each magnetic field line, i.e. 
M(*S>) is a function of two transverse coordinates. 

In particular, starting from an axially (helically) symmetric isotropic plasma equilibrium (|2,1[) . 
one may obtain a family of axially (helically) symmetric anisotropic plasma equilibria (|2.3p , where in 
(|4.ip ^ is the flux function solving the GS (JFKO) equation. Moreover, if the given isotropic plasma 
equilibrium has magneric field lines that are closed loops or go to infinity, one may also obtain a 
wider class of anisotropic plasma equilibria (|2.3p . by choosing the value of the arbitrary functions 
M(^) in (|4.ip separately on every field line. In this case, the resulting anisotropic equilibria may 
have no geometrical symmetries (symmetry breaking). 



4.3 Example of an exact solution: an anisotropic plasma vortex. 

In [10], Bobnev derived a localized vortex-like solution to the static isotropic MHD equilibrium 
system f)2. 1 1) in 3D space. The solution is presented in spherical coordinates (p, 9, (/>) and is axially 
symmetric, i.e. independent of the polar variable 4>. 
The solution has the form 



B = e p V(p) cos 9 + e e W(p) sin 9 + e^U (p) sin 9, P = P - p{p) sin 2 9, (4.6) 

1 Bobnev's solution was found by an ad hoc method without using the GS equation. Since the solution is axially 
symmetric, it corresponds to some flux function ^(r, z) that satisfies the GS equation. However, the explicit form of 
such flux function is not known. 
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where 



U(p)=X n pV(p), p(p)=7P 2 V( P ), W{p) 
V {2X n p) - V (2X n R) 



-V(p) - pV'(p)/2, 



V(p) = B - 



V (x) = 3 



1 - V (2X n R) 



sin x cos x 



7 = Bq- 



V (2\ n R) 



const, Bq,Pq = const. 



l-V (2X n R) 

Here A n ,, n = 1, 2, ... is any member of the countable set of solutions of the equation 

(3 - 4R 2 X 2 n ) sin(2iU„) - 6RX n cos(2i?A„) = 0, R = const. 

(In particular, RXi « 2.882, RX 2 ~ 4.548, i?A 3 » 6.161.) 
Bobnev's solution satisfies boundary conditions 



(4.7) 



(4.8) 



B| p= r — 0, P\ p =r — Po, B 



p=0 



-Boe z , 



(4.9) 



i.e. the domain W S M 3 where the magnetic field is nonzero is a sphere of radius R. The pressure 
outside of W is constant: P = Pq. The magnetic surfaces \& = const inside W are families of 
nested tori of non-circular section, separated by spherical separatrices on which the pressure is 
P = Pq = const; the number and mutual position of the families depends on the choice of R, X n . 

We take R = 1, A = A 3 » 6.161, B = 100, and P = 4500, and find 7 » -72.831. Magnetic 
surfaces P = const and levels of constant magnetic energy density B 2 /2 = const (with P and B 
given by (|4.6p ) are shown in Figure 1 A, B respectively. The spherical separatrix magnetic surfaces 
have approximate radii p\ w 0.376, pi w 0.597. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 



Figure 1: Magnetic surfaces P = const (A) and levels of constant magnetic energy density B 2 /2 = 
const (B) in the static isotropic magnetic vortex (|4.6p (R = 1, A = A3 w 6.161, So = 100, and 
Po = 4500). Here z is the symmetry axis, and x is the radial axis. 
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We apply transformations (|4,ip to this isotropic (MHD) solution, using the arbitrary function 

M(*) = 1 + — sin I — j , #1, * 2 = const, (4.10) 

which is constant on the magnetic surfaces of the original static MHD configuration given by (|4.6p . 
The arbitrary constants *$>i, ^2 are chosen so that M(^) is separated from zero. [In the particular 
example below, we take ^1 = 200, ^2 = 60.] As the result, we obtain an explicit anisotropic (CGL) 
plasma equilibrium configuration {B',p^_,pj|} given by (|4.ip . with the same set of magnetic surfaces. 
The anisotropic pressure components p' ± ,p'^ are no longer constant on these surfaces. Contour plots 
of p' ± ,p'n, and their profile along the radius of the vortex in the direction perpendicular to z, are 
shown in Figure 2. 




Figure 2: Surfaces of constant level of anisotropic plasma pressure components pn (A) and p± (B) in 
an anisotropic vortex (i? = 1, A = A3 ~ 6.161, Bq = 100, and Pq = 4500). Here z is the symmetry 
axis, and x is the radial axis. 

Figure 3(A) shows profiles of magnetic energy densities E = B 2 /Bq and E' = B' 2 /Bq of the 
original isotropic and the new anisotropic plasma vortex respectively. Figure 3 (B) shows profiles 
of normalized anisotropic pressure components p±/Bq ,p\\/Bq along the radius of the spherical 
vortex in the direction perpendicular to the axis of symmetry z. 

The presented solution is an explicit example of a physically meaningful axisymmetric static 
anisotropic plasma equilibrium configuration in 3D space, arising from the model (|2.3p . (|2.5p . This 
solution is regular in the whole domain (a ball of unit radius) and satisfies boundary conditions 

B'| p= i = 0, p' ± \ p =i = p\\\ P =i = const. 

The solution describes a static anisotropic plasma vortex confined by external gas pressure. Other 
choices of the arbitrary function M(\P) lead to different resulting anisotropic configurations. 
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Figure 3: (A) Comparison of normalized energy densities E = B 2 /i3 2 (thick) and E' = (B') 2 /.B 2 
(thin) for isotropic and anisotropic plasma vortex. (B) Comparison of normalized plasma pressure 
components P\\/Bq (thick) and p±/B^ (thin) in an anisotropic vortex in the radial direction. The 
axis x is chosen in the radial direction perpendicular to symmetry axis z. Here R = 1,A = A3 ~ 
6.161, B = 100, and P = 4500. 



5 Conservation law analysis of static MHD and CGL equations 

5.1 Construction and interpretation of conservation laws 

A local conservation law is a continuity equation 

<9 t $° + div$ = (5.1) 

holding locally in 3D space for some physical variable. For dynamical plasma equations, one basic 
example is conservation of mass where <I> = p is mass density and = pv is momentum vector 
density in terms of the fluid velocity v. When static equilibria are considered, such conserved 
densities <3?° are manifestly time-independent, dt&° = 0, while the associated flux vectors 3> are 
divergence free, div$ = 0. In general, a local conservation law of static MHD systems will be a 
vector density <1? which depends on the spatial coordinates x, pressure P, magnetic field B, and their 
partial derivatives with respect to x, such that it is divergence free for all static equilibria. Vector 
densities are physically trivial if they identically have the form of curls, = curl 0, when evaluated 
on static equilibria, with being a local function of the same variables as <1>. A conservation law 
on static equilibria is thus nontrivial if is not such a local curl, namely div $ vanishes essentially 
as a consequence of the static field equations. Conservation laws that differ by a trivial one are 
considered to be physically equivalent. 

The integral of a static conservation law in any domain V in 3D space physically describes the net 
flux through the boundary dV. In particular, if V is a connected 3D region enclosed by a smooth 
closed surface S, then the net flux 



(5.2) 



through S vanishes on smooth static equilibria, while non-vanishing net flux for a static conservation 
law would indicate the presence of a singularity in the flux inside V. More generally, this net flux 
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is independent of S, Si) = J 7 ^, S2), due to Gauss' divergence theorem <f s &-dS—§g 2 <&-dS = 
jy div&dV = where V is the region bounded by the two surfaces Si, S2 in 3D space. 

There is a computational algorithm (Direct Construction Method) for finding conservation laws 
(|5.ip [13,19,31]. For static MHD systems, it is outlined as follows. Because static MHD systems 
are Cauchy-Kovalevskaya type PDE systems (i.e. each system can be written in solved form with 
respect to a first partial derivative of any one spatial coordinate), all of their nontrivial conservation 
laws arise from multipliers whose summed product with the static field equations is identically a 
divergence. Specifically, for isotropic models without fluid flow as considered hereafter, multipliers 
r,A = (A 1 , A 2 , A 3 ) are functions of the spatial coordinates x, pressure P, magnetic field B, and 
their partial derivatives with respect to x, such that 

T(div B) + A • (curl B x B - grad P) = div $ (5.3) 

holds identically (i.e. off of solutions). 

Determining equations for multipliers can be obtained from the fact that divergences (|5.3[) are 
annihilated by variational derivatives (Euler operators) with respect to B and P. In the simplest 
case of multipliers with no dependence on partial derivatives of B and P, the determining equations 
are equivalent simply to the adjoint of the symmetry determining equations f|3.3|) : 

div A = 0, grad T = A x J - curl(A x B) (5.4) 

holding for all static equilibria. For T and A depending on partial derivatives of B and P, there are 
additional determining equations involving variational derivatives of the multipliers themselves. 
The complete system of multiplier determining equations can be solved by an analog of Lie's 
algorithm for solving the symmetry determining equations. When a set of multipliers is known, 
several approaches can be used to compute the flux vector (for a comparison, see [30]). In 
particular, there exists an integral formula involving a homotopy scaling of the field variables [18,19]. 
Alternatively, since the static MHD system (|2.ip admits a scaling symmetry <|3.6j) a. the integral 
formula for <I> can be replaced by a purely algebraic expression [31] in terms of T and A. 

5.2 Conservation laws of the static isotropic plasma equilibrium system 

For static MHD equilibria (|2.ip . we will determine nontrivial conservation laws. In particular, we 
seek multipliers multipliers T and A, such that the summed product of these multipliers with the 
static MHD equilibrium equations (|2.ip yields a nontrivial divergence 

T div B + A • (curl B x B - grad P) = div * (5.5) 

which vanishes when evaluated on static equilibria. An application of the Direct Construction 
Method [19] leads to the following results. 

Theorem 5. The complete set of conservation laws admitted by the static isotropic plasma equilib- 
rium system (|2.ip . for multipliers T, A linear in first partial derivatives of (B, P) and with otherwise 
arbitrary dependence on (x,B,P), is given in TableUl 

Table 1: Conservation laws (15. 5p of the static isotropic (MHD) plasma equilibrium system (|2.ip 
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# 


Multipliers 


Conservation law 


Remarks 


1 
1 




AUr{ f T^ n 

CllV^ • JL J — U 


Conservation of stress and angular stress, 
depending on a general Euclidean Killing vec- 
tor £ = a + b x x, with a, b £ R 3 constant vectors. 


2 


r = /(p), a = -/'(p)b 


div(/(P)B) = 


Conservation of magnetic flux depending on an 
arbitrary function that is constant on magnetic 
surfaces, which reflects the fact that P = const 
on magnetic surfaces. 




p — n A — f 'C P"! T 


HiirC ft P\ 1\ — fl 


£}Y"iL}i"olii7£}/~l \£ i T"/"*"h Ti /"^flr ' o /"*n t' r m t~i lour f"i~f~l nr*~i~ o 
Vjrcllcl cLllZ(_Ll J.Y 11 L-IlllUll b CUIItrilL IdW . rvfcillfcJCLo 

the fact that plasma electric current density J is 
tangent to surfaces of constant pressure. 



In Table [H the symmetric conserved tensor T is the sum of electromagnetic and fluid stress tensors: 

T = -B0B + (P + ^|B| 2 )I. (5.6) 

The arbitrary function f(P) in conservation laws #2 and #3 has, in general, the form / = 
f{fy) given by an arbitrary function that is constant on magnetic field lines (and on magnetic 
surfaces, when they exist). Conservation law #3 involving plasma electric current is analogous to 
conservation of vorticity in time-independent incompressible Euler equations of fluid motion (|2.2p . 

5.3 Conservation laws of the static anisotropic plasma equilibrium system 

For static anisotropic (CGL) equilibria (|2.3p . (|2.5p . we likewise determine multipliers II, ft, and 
T, such that their summed product with the corresponding equations (|2.3p . (|2.5p is a nontrivial 
divergence: 

IldivB + ft ■ ^(1 - t) curlB x B - grad p x - ^rgrad|B| 2 ^j + T(B • gradr) = div*. (5.7) 

Such divergence expressions vanish on anisotropic static equilibria and thus yield conservation laws. 
The following theorem is obtained by an application of the Direct Construction Method [19]. 

Theorem 6. The complete set of conservation laws (15, 7h admitted by the static anisotropic plasma 
equilibrium system (|2.3p . (|2,5p . for multipliers II, ft, T linear in first partial derivatives of (B,p±,r) 
and with otherwise arbitrary dependence on (x, B,pj_,r), is given in Tabled 



Table 2: Conservation laws (15. 7p of the static anisotropic (CGL) system (|2.3p . (I2.5P 



# 


Multipliers 


Conservation law 


Remarks 


1 


n=(l-r)(B-C), 

ft = c, 

T = — B ■ C 


div(C • S) = 0, 


Conservation of stress and angular stress, 
depending on a general Euclidean Killing 
vector C, — a + b x x, with a, b <E R 3 
constant vectors. 


2 


n = /, n = -/ p B, 

t = u + y P \B\ 2 , 

f = f(p,r) 


div(/(p,T)B) = 


Conservation of magnetic flux depending 
on an arbitrary function that is constant 
on magnetic surfaces, which reflects the 
fact p,r = const on magnetic surfaces. 


3 


n = o, n = /'(p)A, 

T = -i/'(p)A.B; 
here A = curl — rB 


div(/(p)A) = 


Conservation of flux related to vorticity of 
the vector field s/1 — rB. 
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In Table O similarly to the isotropic case, the symmetric conserved tensor § is a sum of the 
electromagnetic stress tensor and the anisotropic fluid stress tensor: 

S = —(1 — t)B ®B + (p + i(l — r)|B| 2 )L (5.8) 

The arbitrary function / in conservation laws #2 and #3 is, in general, equal to a constant on 
magnetic field lines (and on magnetic surfaces, when they exist). 

The classification of multipliers and fluxes admitted by the CGL equilibrium system in Table [2] 
coincides with that of the isotropic MHD system in Table [H with the vector field \f \ — rB corre- 
sponding to the isotropic equilibrium magnetic field B and with p corresponding to the isotropic 
pressure P. Note that this relation directly manifests the equivalence of the systems stated in 
Theorem [TJ 



5.4 Example 1: conservation laws of an isotropic plasma vortex (14. 6 j) 



As a simple example, we consider integral forms of conservation laws for smooth static isotropic 
plasma equilibria (Table [1]). If V is a connected 3D region with a smooth boundary S, then integrals 
of conservation laws #2 and #3 in Table [U are respectively 



<£ C • T • n dS = 0, <f f(P)B -ndS = 0, <f f(P)J • n dS = 0, 
Js Js Js 



(5.9) 



IS JS 

where n is an outer normal to S. 

We now write down flux expressions for the Bobnev's isotropic plasma vortex solution presented 
in Section B~3l when S is a sphere of radius p < R = 1. In spherical and cartesian coordinates, the 
arbitrary constant vectors a, b in a Euclidean Killing vector £ = a + b x x have the form 



a = a p e p + age e + 



Qx^x ~t~ ^y^y ~t~ ^e 



flp = «£• sin 6 cos 4> + a y sin 8 sin 4> + a z cos ( 
a e = a x cos 9 cos + a v cos 6 sin <pa z sin 6, 



-a x sin 4> + a y cos 



(5.10) 



with similar expressions holding for b. [Note that a x ,a y ,a z and cartesian components of b are 
constants, whereas spherical components of a and b depend on spherical angles.] 
Using the solution (|4.6p . we find flux expressions 



C-T-n 



-P + p(p) sin 2 6 + - {V 2 {p) cos 2 6 - {U 2 {p) + W 2 {p)) sin 2 9) 



+V(p) sinO cos 6 [U(p)(a^ - rb g ) + W(p){a B +rb$)] 
f(P)B-n = f(P)V(p) cos 6; 

/(P)J-n = 2/(P)^cot0. 

P 



(5.11) 
(5.12) 

(5.13) 



In particular, on each of the two separatrix spheres with radii p = pi,P2 (see Section l4.3p . one 
has V(p) = 0, therefore U(p) = p(p) = 0, W(p) = —pV'(p)/2, and P = Pq = const. Hence on 
these separatrix spheres, fluxes (I5.12p . (|5.13p of the conservation laws #2 and #3 in Table [U vanish 
identically, and the flux (|5.1ip of the conservation law #1 becomes 



C-T-n 



-Po + ^Wsin 2 fl 



P = Pl,P2- 



Substituting (|5.10p . it is easy to verify that the integral f Q * Jq(C ■ T • n) p 2 sin9 d6d(f) vanishes 
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5.5 Example 2: conservation laws for axially symmetric plasma equilibria 

As a second example, we consider conservation laws admitted by static isotropic (MHD) equilib- 
rium system f)2. 1 1) (Tabled]), and explicitly compute fluxes of these conservation laws for axially 
symmetric (Grad-Shafranov) plasma equilibria (|4.3|) satisfying (|4.2|) . 

Since the position vector in cylindrical coordinates (r, (j>, z) is x = re r + ze z , the Euclidean Killing 
vector takes the form £ = a + bxx=(a r + zb^)e T + (a^ + b z r — b r z)e c j ) + (a z — 6^r)e z . Here 

a — d r G r -j- dfjj&fp -\- d z Q z — d x G x ~\~ (XyGy ~\~ CL Z G Z ^ 

b = b r e r + b^e^ + b z e z = b x e x + b y e y + b z e z , 

are arbitrary constant vectors in M 3 ; 

a r = a x cos (j) + a y sin <fi, = a y cos <fi — a x sin <fi, 
b r = b x cos 4> + b y sin (j), b^ = b y cos 4> — b x sin <^>; 

where a x ,y,z, t>x,y,z = const. 
Denoting the total plasma energy density e = P(*&) + ^|B| 2 , and magnetic field components 

^ z iNf) * r 
B = 5 r e r + Bses + 5 2 e 2 = — e r H e 2 , 

we find the following six conservation laws originating from conservation of stress and angular stress 
in Table CD 

1 d 1 d 

-q^ [r(e - B 2 ) cos <j> + rBj-B^ sin 4>] - - — [(e - -B 2 ) sin + B r B^, cos 0] 

- |- [B z (B r cos - ^ sin 0)] = 0, (5.14) 

oz 

Id Id 

-7T k(e - -Br) sm <^ ~ TB r B < h cos 0] + -— [(e - .B|) cos(/> - -B r -B<A sin 0] 
r Or 7™ acp J 

-J^£*(J3 r sin0 + ^cos0)]=O, (5.15) 

^ [rBri?J -^[ e -^ =0 ' (5 - 16) 
13 1 3 

— |r(z(e - B 2 ) + rB r B z ) sin </> - rzB r BA, cos 0] + -— [z(e - £??) cos 4> + (rB z - zB r )Bs sin 
r or r ocp v 

- y z [( r ( e - B z) + zB rB z ) sin + zB^B z cos 0)] = 0, (5.17) 

I Q I Q 

-TT \r(z(e - B 2 ) + rB r B z ) cos <fi + rzB r BA, sin d>] — \z(e - B"i) sm<p - (rB z - zB r )B^ cos 

r or r ocp v 

- ^[{r(e - B 2 Z ) + zB r B z ) cos - zB^B.smcJ))} =0, (5.18) 
1 d d 

[r 2 B r B^\ + — [rB^Bz] = 0. (5.19) 
r or Oz 

For axially symmetric plasma configurations, the conservation laws corresponding to conservation 
of magnetic flux and vorticity in Table Q] become trivial. 

Remark 2. Conservation laws (I5.16h and (I5.19h have the axially invariant form J^<I>i(r, z) + 
J^&2( r , z) = and are admitted by the Grad-Shafranov equation (|4,2p for any choice of arbitrary 
functions I(^f), P(^f). The other four conservation laws (|5,14p . (|5.15p . (|5.17p . (|5.18p are admitted 
by the full static MHD equilibrium system (|2.ip but not by the Grad-Shafranov equation (|4.2|) . 
since they explicitly contain the angular variable (p. 
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6 Conclusions 



Magnetohydrodynamics (MHD) and Chew-Goldberger-Low (CGL) models are the two most widely 
used continuum plasma descriptions, valid for the cases of isotropic and strongly magnetized 
(anisotropic) plasmas respectively. Knowledge of analytical properties of these nonlinear PDE 
systems and, in particular, methods of finding exact solutions, are highly important for applica- 
tions. 

In this paper we have classified complete sets of admitted point symmetries and conservation laws 
of the static isotropic plasma equilibrium system (|2.ip and the static anisotropic plasma equilibrium 
system (|2.3h . (|2.5p . This classification has led to establishing a direct transformation (|2.6[) . (12. 7p 
between the two systems. The transformation implies the equivalence of solution sets: every 
static anisotropic (CGL) equilibrium can be obtained from a solution of static isotropic (MHD) 
equilibrium system, and to each solution of the static MHD system there corresponds an infinite 
family of CGL equilibria depending on an arbitrary function defined on the set of magnetic surfaces. 

The established equivalence yields an effective procedure of construction of exact explicit anisotropic 
(CGL) static plasma configurations from a single known MHD or CGL solution. Many physically 
meaningful static MHD solutions are known, and each of them gives rise to families of anisotropic 
(CGL) equilibria with the same topology of the magnetic field. It follows that all axially and he- 
lically symmetric CGL configurations can be found from solutions of conventional Grad-Shafranov 
and JFKO equations. An example of an explicit solution describing an anisotropic axially symmet- 
ric plasma vortex is given in Section HI 

We note that symmetry classification of the static isotropic MHD system ()2. 1|) does not lead to a 
symmetry classification of the Grad-Shafranov or JFKO equations, since the latter are reductions 
of the system (|2.ip . The conservation laws found in this paper for the static isotropic MHD system 
(12. 1|) yield some particular conservation laws of the Grad-Shafranov equation (Section 15. 5h . but 
again not a full classification. 

A natural next step will be to accomplish a similar complete symmetry and conservation law 
analysis of dynamic (v ^ 0) plasma equilibrium models, and further, of time-dependent (non- 
equilibrium) equations. Dynamic equilibrium models are already known to possess rich symmetry 
structure [20,32], but the complete analysis was not done due to the complexity of these PDE 
systems. Recently developed symbolic computation software [21, 22] used in this paper will be 
applied to study symmetries and conservation laws of these systems. 
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